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Abstract 

A linear relation for charged-particle emissions is presented starting from the microscopic mech- 
anism of the radioactive decay. It relates the logarithms of the decay half-lives with two variables, 
called x' p\ which depend upon the Q-values of the outgoing clusters as well as the masses and 
charges of the nuclei involved in the decay. This relation explains well all known cluster decays. 
It is found to be a generalization of the Geiger-Nuttall law in a radioactivity and therefore we 
call it the universal decay law. Predictions on the most likely emissions of various clusters are 
presented by applying the law over the whole nuclear chart. It is seen that the decays of heavier 
clusters with non-equal proton and neutron numbers are mostly located in the trans-lead region. 
The emissions of clusters with equal protons and neutrons, like ^^C and ^^O, are possible in some 
neutron-deficient nuclei with Z > 54. 

PACS numbers: 21.10.Tg, 23.60. +e, 23.70.+j, 27.60.+j, 27.90. +b 
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I. INTRODUCTION 



Charged-particle emissions are among the most important decay modes of atomic nuclei. 
Almost all observed proton-rich exotic nuclei starting from A ~ 150 are a radioactive P]. A 
substantial number of proton decays have been observed in proton-drip-line nuclei around 
the rare earth region 2|. The spontaneous emission of charged fragments heavier than the a 
particle (cluster decay) was predicted in Ref. [3] and later established experimentally in trans- 



flag. 



lead mother nuclei decaying into daughters around the doubly magic nucleus ^ospi^ 
Even a second island of cluster radioactivity was predicted in trans-tin nuclei decaying into 
daughters close to the doubly magic nucleus ^°°Sn [sl. 



A number of theoretical models were prop osed to describe t 
|9|,[lo|,[ll|,[l2l,[l3 



process 



3, 15, 3, 17, 3, 3 (see also Refs. 



le chargec 
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-particle decay 



24 



25| for very 



recent calculations). In general the decay process, ranging from proton to heavier cluster 
radioactive decays, can be described by a two-step mechanism [2^. The first step refers to the 
formation of the particle and its motion on the daughter nuclear surface. In the second step 
the cluster, with the formation amplitude and corresponding wave function thus determined, 

is assumed to penetrate through the centrifugal and Coulomb barriers l27l . l28l . y,0. This 



second step is well understood since the pioneering work of Gamow 27|. In macroscopic 



models, cluster deca y is treated as_ the quantum tunneling process of an already preformed 
particle 
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25[ |. where features like the probability that the cluster is 



formed on the nuclear surface are ignored. In these models the clusterization process is 
included in an effective fashion by introducing quantities adjusted to reproduce as many 
measured half-lives as possible. Such semiclassical models are being successfully applied 
even at present, although in some cases microscopical ingredients are also included IJ]. In 
microscopic theories the formation amplitude is evaluated starting from the single-particle 
degrees of freedom of the neutrons and protons that eventually become the cluster. This is 



generally a formidable task which requires advanced computing 



theoretical schemes to describe the clusterization process 15 
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aci 
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ities as well as suitable 



18|. 



On the other hand, this variety of theoretical models may serve as a guide to our searching 
of semiclassical relations in the radioactive decay. The first striking correlation in a decay 



systematics was noted by Geiger and Nuttall 



3l| . This relates the decay half-lives T1/2 and 
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decay energies Qa as, 

logTi/2 = aQ;i/2^6, (1) 

where a and h are constants. Nowadays it is understood that the Q-value dependence in 
Eq. ([1]) is a manifestation of the quantum penetration of the a-cluster through the Coulomb 
barrier (see, for example, Ref. [32|]). But this equation ignores the probability that the a- 
particle is formed on the nuclear surface starting from its four constituent nucleons moving 
inside the mother nucleus. The linear relation ([T]) has been found to hold well for the ground- 
state to ground-state decays of even-even nuclei in the same major shell with fixed proton 
number. However, the Geiger-Nuttall law in the form of Eq. ([1]) has limited prediction power 



since the coefficients a and h change for the decays of different isotopic series 



331]. Intensive 



works have been done trying to generalize the Geiger-Nuttall law for a universal description 



of all observed a decay events [33 
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381]. For example, in the work of Viola and 



Seaborg [33|, the a and h coefficients of Eq. ([T]) are assumed to be linearly dependent upon 
the charge number of the daughter nucleus. But the physical origin of this dependence is not 



clear. Empirical linear relations were also found in the cases of proton decays 
heavier cluster decays 2^, 
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441 ] and both a and heavier cluster decays 



3£ 
45 



46 



37 



4q|,^], 
47]. 



461 ] and thus will not 



Reviews on existing empirical relations can be found in Refs. 
be detailed here. In particular, some recent searches of correlations in radioactive decays 
start from the macroscopic description of the decay process with a concentration on the 



Coulomb barrier penetrability 



32 



34 



37 



45| . which are physically more sound than mere 
empirical relations. But in these macroscopic approaches one has to assume an effective 
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Besides, an effective spectroscopic 



interaction between the cluster and the core 
factor (formation amplitude) has to be introduced 37|] which, however, is model-dependent 
and sensitive to details of the effective interaction. 

In a recent Letter 48|, we have introduced a linear universal decay law (UDL) starting 



from the microscopic mechanism of the charged-particle emission. Our aim was to find a 
general framework valid for all clusters, which may be used in the future as a gauge to probe 
effective formulas. This is an interesting subject in itself, but perhaps even more important 
is that it may help in the ongoing searching of new cluster decay modes from superheavy 
nuclei 29|] . The UDL relates the half-lives of monopole radioactive decays with the Q- values 
of the outgoing particles as well as the masses and charges of the nuclei involved in the 
decay, reflecting quite well the systematical trend of experimental data. In this paper we 



will complete the brief presentation given in Ref. 48| with details of the construction of the 
formula and approximations leading to it. Besides, we present the predictions of the UDL 
on the most likely emissions of various clusters. 

In Section II is the Formalism. In Section III systematics of experimental a and cluster 
decay half-lives are analyzed and compared with the corresponding calculations. In Sec- 
tion IV possible observations of new cluster decays are suggested. A Summary and the 
Conclusions are in Section V. 



II. FORMALISM 



In a classic paper 26|], Thomas derived the expression of the cluster decay width by 
evaluating the residues of the corresponding S-matrix in the framework of the R-matrix 
theory 49|]. The decay half-life thus obtained has the form, 

2 



7l/2 



^ln2 In 2 



(2) 



where is the decay width, v the outgoing velocity of the charged particle carrying an 
angular momentum /. R is the distance between the corresponding centers of mass of the 
cluster and daughter nucleus, which should be large enough that the nuclear interaction 
is negligible. is the Coulomb- Hankel function and its arguments are standard, i.e., 

p = jjiuR/Ti and the Coulomb parameter is x = 2ZcZd(? j'hv with p being the reduced mass 
of the cluster-daughter system and and Z^ the charge numbers of the cluster and daughter 
nucleus, respectively. Eq. ([2]) contains the two-step mechanism mentioned above. The 
quantity Fc{R) is the formation amplitude of the decaying particle at distance R, which is 
usually evaluated as the overlap between the mother wave function and the antisymmetrized 
tensor product of the daughter and cluster wave functions. The penetrability is proportional 
to \Hj^ {x, p)\~'^ ■ This equation is the basis of all microscopical calculations of radioactive 



decay processes jiol . l39l ]. It is valid for all clusters and for spherical as well as deformed 



cases. The ratio Ni = RFc{R) / (R) , and therefore the half-life itself, is independent of 



the radius R 



16|. In Ref. [5o| it is shown that the expression of Eq. ([2]) coincide with the 



quantum-mechanical interpretation of the half-life as the outgoing flux per unit of time. 

In what follows we will apply the exact expression of Eq. ([2]) . Our aim is to find a few 
quantities that determine the half-life. Expanding in these quantities we hope to be able to 
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find, at the lowest order of perturbation, an expression of the half-life which is as simple as 
the Geiger-Nuttall law but valid in general, i.e., for all isotopic series as well as all type of 
clusters. This is possible since Eq. itself is valid in general. The number of variables that 
we have to look for should be small for cases of interest, i.e., for the decay of medium and 
heavier nuclei. In fact most interesting is the predicting power with respect to superheavy 
nuclei, which are at the center of attention of present experimental activities. With this in 
mind we notice that the Coulomb- Hankel function can be well approximated by an analytic 
formula, which for the / 



5ii 



channel reads 
Ho{x,p) ~ {cotpy/^exp [x(/5-sin/3cos 
where the cluster Q- value is Qc = and 



2 o P 

cos p = — 

X 



QcR 



(3) 



(4) 



One sees that cos^/3 would be a small quantity if Z^Z^ is large. In this case one can 
expand the last term in a power series of cos (3 (with (3 = arccos(cos/5)) as, 

(5) 



n ■ ^ ^ ^ n COS^/? COS^/? 



and for medium and heavier nuclei (the heavier the better) terms beyond the third order 
can be neglected. One obtains. 



log 



\H+{x,p)\', 2x 



cot P 



In 10 



1/2 
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3/2 



(6) 



and therefore, 



logTi/2 



2X 
In 10 

+ log 



TC 



1/2 
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cot /3 In 2 
uR^\F,{R)\^ 



+ 



3 (x. 



3/2' 



(7) 



which is dominated by the first two terms. For the radius R in above equation, one can take 
the standard value of = Ro{A]/^ + A]/^) with Ro ~ 1.2 fm |39|, lisj. Defining the factors 
x' and p' as, 



A 



X 

p' 



-X — ZcZfi\l — , 
=2, /2m V Qc 



^/2mRoe^ 



(px) 



1/2 



(8) 



where A = AdAc/{Ad + Ac) = fi/m and m is the nucleon mass (within the errors of our 
treatment we take mc^ ^ 938.9 MeV and he = 197.3 MeV fm), one gets 

, ^ V2Me^7i , 4eV2MR~o , 

= luiuT^ - MnlO P 
, / cot/51n2 \ 

= ax' + 6p'-21og|Fc(i?)| + c, (9) 

where a and h are constants and o(3) corresponds to the remaining small terms in the 
Coulomb penetration. The terms o(3) and log cot /3/(z/i?^) change rather smoothly for the 
decay cases of interest and may be safely approximated as a constant c. 

A straightforward conclusion from Eq. ([9]) is that logTi/2 depends linearly upon x' and 
p' . Still the strong dependence of the formation probability upon the cluster size has to 
be taken into account by Eq. ([9]). This seems to be a difficult task, since the formation 
probability is strongly dependent upon the nuclear structure of the nuclei to be analyzed. 
In other words, if such a simple linear relation is correct one has to be able to demonstrate 
that the formation amplitude depends only linearly upon x', p' or an additional variable. 
We found that this is indeed the case by exploiting the property that for a given cluster 
A^o = RFc{,R) I Hq {xi p) does not depend upon R. Using the approximations leading to Eq. 
([9]) one readily obtains the relation, 

Aey/2M 



log 



R!F,{R') 



RFc{R) 



/ilnlO 



Ro p', (10) 



where R' = R'^iA]^^ + Al^^) is a value of the radius that differs from R. The above equation 
can also be written as, 

log \RF,iR)\ ^ log \R'F,{R')\ + - p'. (11) 

Since for a given cluster any nuclear structure would be carried by the terms RFc{R) and 
R'Fc{R') in exactly the same fashion, Eq. (ITT]) implies that the formation amplitude is indeed 
linearly dependent upon p' . Therefore one can write |48l |. 



logTi/2 = ax' + &p' + c. (12) 

The coefficients h and c in this relation are different from that of Eq. ([9]) since the terms 
hp' + c have to include the effects that induce the clusterization in the mother nucleus. 



This relation holds for the monopole radioactive decays of all clusters and we called it the 



UDL 48|. The relation can be easily generalized to include the / 7^ decay cases by taking 
the effects of the centrifugal potential on the barrier penetrability into account [51^. 

It can be easily recognized that the UDL includes the Geiger-Nuttall law as a special 
case since p' remains constant for a given a-decay chain and x' oc Q'^l"^ . Besides, one basic 
assumption behind the relation (fT2|l is that one can define a proper radius Bl that leading 
to a stable formation amplitude Fc{R') for all cluster radioactivities. In the next Section we 
will probe these conclusions, and the approximations leading to them. 

III. SYSTEMATICS OF EXPERIMENTAL DATA 

In this Section we will analyze ground-state to ground-state radioactive decays of even- 
even nuclei. We take all a decay events from emitters with 78 < Z < 108 for which 
experimental data are available to us. We take the data from the latest compilations of 







Refs. Ijj, n2l and the lists of Refs. 



22 



531 1 . For the decay of heavier clusters we have selected 
11 measured events ranging from ^^C to ^'^Si [7|. The branching ratios of the cluster decays 
relative to the corresponding a decay are in the range of 10~^ to 10^^^. The partial half-lives 
of observed cluster decays are between lO^^s and lO^^s. 

A. Experimental constraint on the formation amplitude 

The formation amplitude Fc{R) reflects the nuclear structure effect on the cluster de- 
cay process. According to Eq. ([2]) the formation amplitude Fc{R) can be extracted from 
experimental data as, 

log \RF^{R)\ = -\ logTf/f • + 1 log ■ 



V 



■K(X,P)P 



(13) 



By using Rq = 1.2 fm we evaluated the function log \RFc{R)\ corresponding to a clusters 
to obtain the results plotted in Fig. [H One sees that as a function of the charge number 
of the emitters the formation probabilities are located in the range log \RFc{R)\ = —1.5 ~ 
—0.75 fm~^/^ with about half of the data below -1.0 fm^^/^. The formation amplitude is 
therefore in the range Fc{R) = (0.03 ~ 0.18)/-R (where R is in fm and Fc{R) in fm~^/^). 
One thus confirms that for a given cluster the formation amplitude is constant within an 
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FIG. 1: The a decay formation amplitudes log \RFc{R)\ as a function of the charge number of the 
mother nucleus Z . 

order of magnitude. The stability of the a decay formation amplitude indicates that the 
linear relation described by Eq. ( |T2l) is not as unexpected as one might have assumed. 

It is seen from Fig. [1] that in a few cases the formation amplitudes become small, i.e., 
with the log I i?-Fc(-R) I < —1.4 fm~^/^. These correspond to the a decays of nuclei ^^^Pb, 

error in the experimental half-life of this nucleus is 
still large . The a formation amplitudes may have been significantly reduced in the 

former four nuclei which are approaching the Z = 82 and/or = 126 shell closures. 

Following the same procedure as above we evaluated log |i?Fc(i?)| for observed heavy 
clusters, as seen in Fig. [2l One sees that now log |i?Fc(-R)| is in the range -9 to -3, i.e., 
FjyK) = (10~^ — 10^^)/_R fm"^/^. Given the variety of clusters (from ^"^C to ^"^Si) involved 
in the Figure this wide range of 6 orders of magnitude is also expected. 

We are now in a position to probe the validity of the linear relation between the logarithm 
of the formation probability as a function of p' as implied by Eq. fill I) . As seen in Fig. [3] 
that relation holds rather well. The majority of available experimental data corresponds to 
a-decay. Since the formation amplitude of a given type of cluster is rather constant one 
sees in the Figure an accumulation of black points around log \RF(R) \ ^ —1.2 fm^^/^. This 
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FIG. 2: Same as Fig. [T]but for those of heavier cluster decays. 
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FIG. 3: (Color online) The formation amplitudes log\RFc{R)\ for both a and cluster decays as a 
function of p' . 



9 



corresponds to a-decay events. For the other clusters the formation probabihties decrease, 
as expected. The important point that one can make from the Figure is that the predicted 
hnear trend is confirmed. The hnear trend of experimental data ensures that one can find 
an radius R' where Fc{R') remain constant for all radioactivities. The R' value can be 
determined through a fitting procedure. 

There is a deviation to this trend at p' ~ 460 which corresponds to the decay of 
^^'^U — > 206jjg _j_ 28]y[g_ may expect that the formation of ^^Ne is more favored in 

^^^U (with the daughter system of ^°®Pb). However, the decay emitting ^^Ne is hindered due 
to the presence of a much lower Qc value (i.e., larger %'). 



B. Systematics with the UDL 

The prediction power of the UDL (Eq. ( |T2i) ) on radioactive decay of medium and heavier 
nuclei has already been shown in our previous Letter [48]. Essentially, only the coefficients 
b and c are free parameters not provided by the UDL. Without loss of generalization, in 
that paper all coefficients of Eq. f|T2|) are determined by fitting experimental data. The 
inclusion of a as a free parameter takes into account the effect of higher-order terms of the 
Coulomb penetrability. For example, the constants a, b and c corresponding to a-decay were 
determined to be 0.4065, -0.4311 and -20.7889, respectively [48i]. The standard root mean 
square (rms) deviation between the UDL and experimental a-decay half-lives is cr = 0.3436. 



For the rms deviation we take the definition of Ref. 



37| 



n 



— E[iog(^''^V^'^'^^")]' , (14) 



where n is the number of decay events included in the fit and T^^p*' and T^*^'' the ex- 
perimental and calculated decay half-lives, respectively. The fitted value for the coef- 
ficient a is close to the value calculated by its definition in Eqs. (Q and f|T2l) . namely 
a = e^7rV2m/ {hln 10) = 0.4314. Even if we fix the coefficient a at this value, the descrip- 
tion power of the UDL is still encouraging, as illustrated in the upper plot of Fig. HI With 
only two free parameters of b and c, the UDL can reproduce experimental a decay half-lives 
with a rms deviation of cr = 0.4606. 

In some cases the experimental half-lives are noticeably underestimated by the UDL, with 
yExptyyCai. > 4_ rj.^^^^ correspond to the a decays of ^^^'^lopb, 208,2iopQ ^^^^ 2i2j^^_ rj.^:^^ 
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TABLE I: Upper: Coefficient sets of Eq. (jl2|) that determined by fitting to 



experiments of a decays 



481 ] ■ and the corresponding 



(I), cluster decays (II) and both a and cluster decays (III), respectively 
rms deviations; Lower: same as the upper part but with coefficient a fixed to its calculated value 
of a = 0.4314. 







Il(cluster) 


III(a+cluster) 


a 


0.4065 


0.3671 


0.3949 


b 


-0.4311 


-0.3296 


-0.3693 


c 


-20.7889 


-26.2681 


-23.7615 


a 


0.3436 


0.6080 


0.6107 


a 


0.4314 


0.4314 


0.4314 


b 


-0.4608 


-0.3921 


-0.4087 


c 


-21.9453 


-32.7044 


-25.7725 


a 


0.4606 


0.7901 


0.7631 



deviation may be related to the fact that the formation amplitudes in these nuclei, due to 
the shell closures of = 126 and Z = 82, are significantly smaller than those in the open 
shell region. 

In fig. m we plotted calculations with the UDL on a and heavier cluster decay half-lives 
and comparisons with experimental data. In these plots, the coefficient a of the UDL is taken 
as its calculated value while h and c are determined by fitting to corresponding experiments. 
In the figure we plotted the quantity logTi/2 — hp' as a function of x'- Similar linear trend 
can be achieved if we plot experimental data as a function of p', as seen in Fig. [5] where all a 
and heavier cluster decays are considered. For all observed a and heavier cluster decays, the 
x' and p' values are in the wide ranges of 105 < x' < 640 and 60 < p' < 660, respectively. As 
a result, the functions logTi/2 — &p' and logTi/2 — o-x' we plotted in Figs. [Hand [5] change over 
200 orders of magnitude. But the decay half-lives are in the range of —8 < logTi/2 < 28 (in 
seconds) . 

In Table [T] the constants a, h and c that fit the data sets of a as well as cluster decays 
are collected. The fitted values with a as a free and fixed parameter are shown in the upper 
and lower part of the Table, respectively. In the table we also give the corresponding rms 
deviations between experiments and UDL calculations with these coefficient sets. 
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FIG. 4: (Color online) UDL description of a decays (upper panel) and heavier cluster decays 
(lower panel) with the coefficient a fixed at its calculated value. The black points correspond to 
experimental data with decay half-lives given in seconds. The lines are given as ax' + c with c 
values from the lower part of Table [H 

IV. PREDICTIONS AND DISCUSSIONS 

Using the UDL it is straightforward to evaluate the half-lives of all cluster emitters 
throughout the nuclear chart if reliable values of the binding energies (i.e., of the cluster 



Q- values) can be obtained. We do this by using the latest compilation of nuclear masses (52 1 
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FIG. 5: (Color online) UDL description of both a and heavier cluster decays as a function of p' . 
The lines are given as bp' + c with coefficients (set III) from the lower part of Table HI 

With the Q-values thus obtained we evaluated the decay half-hves of all isotopes included 
in that compilation by applying the UDL. We will first show the case of the decay of a 
particles and afterwards that of other relevant clusters. For simplicity in what follows, only 
results calculated with coefficients from the upper part of Table [T] are shown. 

Since the half-lives of decaying nuclei which live a very short or very large time can not 
be measured we will only consider even-even a emitters with half-lives in the 30 orders of 
magnitude range —10 < logTi/2 < 20 (in second). The UDL predictions of the corresponding 
half-lives are shown in Fig. [HI employing the coefficient set I from the upper part of Table [B 
Within the constraints that we imposed a total number of 269 even-even a emitters have 
been found, which have charge numbers Z > 52. It is seen from the Figure that the most 
favored a decays are from neutron-deficient nuclei around the trans-lead and superheavy 
regions. 

Since, as mentioned above, the emitters that we have used to determine the coefficients of 
the UDL (Table [11) have charge number Z > 78, it would be interesting to probe the law for 
nuclei with Z- values below that limit. We thus took the extreme cases of decays from nuclei 
in the trans-tin region. The a decay properties of nuclei in this region has been intensively 
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FIG. 6: (Color online) Prediction of the UDL (with the coefficient set I) on the logarithms of 
half- lives (in seconds), logT^/21 the a decays of even-even nuclei. 



TABLE II: Experimental and UDL calculated values (with the coefficient sets I and III) of a decay 
half-lives (in seconds) of even-even nuclei in the trans-tin region. Experimental data are taken 



from Ref. [l| except for the half-life of ^^'^Xe which is from Ref. 



Emitter 



Qa(MeV) 



logTg^,'-(ni) 



108 rpg 

iiOXe 
ii2xe 



4.290 
3.445 
3.885 
3.330 
3.534 



-4.155 
0.6320 
-0.7850 
2.477 
1.770 



-3.446 
0.9761 
-0.3774 
2.951 
2.861 



-4.484 
-0.1812 

-1.441 
1.799 
1.766 



studied in recent years [55|, l56|, l57| . In Table [TTl we compare experimental results on the a 
decays of Te, Xe and Ba isotopes and the predictions of the UDL. One sees that in all cases 
the experimental values lie between the ones calculated by using the parameters of the sets 
I and III in Table [H confirming the prediction power of the UDL. 

We will now apply the UDL to evaluate the emissions of heavy clusters which are good 
candidates to be observed, namely ^^'^^C, i6,i8,20q^ 20,22,24^^^ 24,26,28^2 and 28,30,32,34gj^ q^, 

served cluster radioactivities exhibit much longer partial half-lives than those of the corre- 
sponding a decays. This can be easily understood if we compare the x' values of the heavier 
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FIG. 7: (Color online) x' values for a and ^^C radioactivities as a function of the mass numbers of 
mother nuclei. 

cluster and a radioactivities since the logarithm of the half-life is proportional to x'- As 
a typical example, in Fig. [7] we plotted the x' values of a and ^^C decays as a function of 
the mass numbers of the mother nuclei. The x' values of heavier clusters are mostly much 
higher than those of the corresponding a decays, indicating that it is more difficult for the 
heavier clusters to penetrate through the Coulomb barrier. Besides, from the figure one sees 
that nuclei favoring cluster decays should mostly be located in the trans-lead region. 

In Fig. [S] we show the predicted half-lives corresponding to the most favored cluster 
radioactivity, namely ^^C decay. Our calculations show that nuclei like ^^°'^^^'^^^Ra, 222,224r]n]^ 
and 226-[j (^g^j^ have partial decay half-lives shorter than 10^^ s, among which the ^^C decays 
of ^^^'^^"^Ra have been observed 0, E]. 

Our calculations also show that nuclei that most probably emit clusters with non-equal 
proton and neutron numbers like ^^C are concentrated in the trans-lead region. This is 
consistent with the expectation from the schematic picture of Fig. [71 For heavier clusters 
the formation probability is even smaller and therefore the corresponding decay probability 
is also smaller. As another typical example, in Fig. [9] we plotted calculations for the half- 
lives of the ^^Ne radioactivity. One sees in this Figure that the shortest half-lives correspond 
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N 

FIG. 8: (Color online) Same as Fig. [6] but for the ^^C cluster decay and with the coefficient set II. 

to mother nuclei around Z=92 and N=138. In all cases this half-hfe is larger than 10^^ s, 
which is many orders of magnitude larger than the cases corresponding to the decay of ^^C 
analyzed above. 

All heavier-cluster-decaying nuclei decay also by emitting a particles. In fact a decay 
is usually the overwhelming dominant decay channel, as seen from Fig. [71 Therefore in 
planning the detection of a probable cluster decay one has to consider carefully the branching 
ratio of the a-decay channel relative to the corresponding cluster decay of interest, i.e., 
^rei. ^ Tf/2/T{=)^'*'"'. The log 6""^' values are negative, which can not be too small for the 
heavier cluster decay to be detectable. We can evaluate these branching ratios by using the 
UDL. To search for probable cluster emitters we select particle decay channels for which 
neither the half-lives are too large nor the branching ratios are too small. We thus use the 
criteria T1/2 < 10^° s and If^^' > 10^^^ which is two orders of magnitude outside present 
experimental limits. The corresponding calculations for the emissions of Nc 7^ Zc clusters 
are listed in Table IIIII To give an insight of the expected precision, we present in Table IIVI 
comparisons between calculations and experiments for the half-lives of the eleven observed 
heavier cluster decay events of Fig. |H 

The a-decay mode dominates the decays of all heavier-cluster emitters we listed in Ta- 
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FIG. 9: (Color online) Same as Fig. [6] but for the Ne decay. 



ble mil In most cases the branching ratio between a decay and all other decay channels 
(including /3 decay) is 6° ~ 100% [1]. But there are exceptions, in particular the nuclei 
^^^Pu, ^^^Pu and ^^'^Cm which have the a-decay branching ratios of 6" = 11%, 6" = 6% and 
6° < 10%, respectively [l]. 

We will now analyze the more rare case of radioactive decays of A^^^ = clusters heavier 
than the a particle. Intense studies have been made in the prediction and searching for the 



emissions of A^^ = Z,. clusters 



56 



57 



58, 



59 



60 



61 



62j. Experiments have not pinned down 



the observation of these clusters yet, although efforts have been made, particularly looking 
for the probable emission of ^^C [60]. We have therefore apply the UDL to investigate 
regions in the nuclear chart where such cluster would be likely to be formed and emitted. 
The half-lives of ^^C decays thus calculated are plotted in Fig. [TOl The emissions of other 
A^c = Zc clusters like ^^O show similar patterns. It is seen from the figure that Nc = Zc 
cluster emitters form two islands, decaying into daughter nuclei around ^°°Sn and ^°^Pb. 
This is consistent with theoretical calculations using the fission model [s], [6l| . 

A first glance at Fig. [10] may suggest that the emissions of Nc = Zc clusters like ^^C 
should be more favored than those of other 7^ Z^. isotopes since the former particle is 
usually more tightly bound. Such a picture is also expected if we compare the x' values for 
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TABLE III: Predictions of the UDL on probable 
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the radioactive decays of other isotopes. A typical example is given in Fig. [TT] where we 
plotted the x' values of the ^^C and ^^C radioactivities. It is seen that the x' values of the 
^^C radioactivity are mostly smaller that those of ^^C, indicating that it should be much 
easier for the ^^C particle to penetrate through the Coulomb barrier, especially in nuclei 
close to the proton drip line. However, the probability of the decay of Nc = Zc clusters 
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TABLE IV: Experimental and UDL calculated values (with the coefficient set II of cluster decay 
half- lives (in seconds). The experimental values are from Ref. Q]- 



Emitter Mode Qc(MeV) logTi/2(II) 



222Ra 14q 3305 11.01 11.07 

224Ra 30.54 15.86 15.59 

226Ra 28.20 21.24 20.33 

228Th 20q 44.72 20.72 21.59 

230u 22]^g gl_39 19 22 20.73 

230Th 24]^g 57_7g 24.61 24.74 

232u 24]sjg g2.31 20.40 20.68 

234u 28]Yig 74.11 25.75 25.36 

236pu 28]Yig 79. Q7 21.52 21.02 

238pu 32g- 91 19 25.27 25.39 

242cm 34si 9g^5i 23.15 22.87 



become small if we take into account the fact that the likely emitters are mostly close to 
the proton drip line and are dominated by the decay mode of With the same selection 
criteria discussed above, our predictions on probable emissions of = Zc clusters are listed 
in Table El Since in all cases the decay by the emission of an a particle is much more likely 
than the corresponding decay by the emission of heavier clusters, in Table |V] we only show 
emitters that are known to decay a It is seen from the Table that the mostly likely 
Nc = Zc cluster emitter is the nucleus ^^^Ba. 

V. SUMMARY AND CONCLUSIONS 

Starting from the exact expression for the half-life of cluster decaying nuclei (Eq. ([2])) we 
found that this expression is dependent upon a quantity called cos^ (3 which for medium and 
heavier nuclei is small (Eq. (jl])). For the case of / = (monopole) transitions we expanded 
the exact expression to the lowest order in cos^ (3 and used the property that the half-life 
does not depend upon the matching radius R (Eq. ffTTl) ). We thus found that the logarithm 
of the half-life is linearly dependent upon two parameters, called x' p' i which depend 
only upon the Q-value of the outgoing cluster and of the charges and masses of the particles 
involved in the decay (Eq. ([8])). The resulting linear expression (Eq. (fT2l) ) is found to be 
a generalization of the Geiger-Nuttall law and we call it the universal decay law (UDL). 
The UDL is valid for all / = transitions. This monopole linear equation contains three 
constants, called a, h and c. We fitted the experimental half-lives of ground-state to ground- 
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FIG. 11: (Color online) x' values for ^^C and ^^C radioactivities as a function of the mass numbers 
of mother nuclei. 
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TABLE V: Predictions of the UDL on probable emissions of Nr 
branching ratios, b°^, of the cluster emitters are taken from Ref. 1| 
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state ct-decay and heavier cluster decay processes in even-even nuclei to obtain the values 
of the constants given in Table [B We found that the UDL predicts with great precision 
the half-lives of radioactive decays, both a- and cluster-decays, and for all isotopic series, 
as expected since the original exact expression for the half-life is valid in general. This law 
may also help in the ongoing search of new cluster decay modes from superheavy nuclei. 

Using the UDL we have evaluated the decay half-lives of various cluster emitters through- 
out the nuclear chart with reliable values of binding energies as input. It is found that the 
a decay is favored in neutron-deficient nuclei around the trans-lead and superheavy regions. 
The decays of heavier clusters with non-equal proton and neutron numbers are mostly lo- 
cated in the trans-lead region. The probability of the decay of clusters with equal numbers 
of protons and neutrons is small since the likely emitters are mostly close to the proton drip 
line and are dominated by the decay mode of (5^ . 

An important conclusion from the UDL is that the cluster formation amplitude Fj^K) 
is exponentially dependent upon the variable p'. The implication of this linear trend on 
nuclear structure effects may deserve further investigations in the future. 
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